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^ Abstract 

In a Hartree-Fock picture, itinerant ferromagnetism results from a competi- 
tion between kinetic and exchange energy, with the magnetized state being 
favored at large interaction strength. In a recent paper Q we showed that in 
O ' contrast to this average effect, fluctuations of off-diagonal interaction matrix 

elements systematically reduce the ground-state magnetization. When the in- 



depends on the ratio A between average exchange energy and the fluctuation 



teraction dominates, the occurence of a non-zero ground-state magnetization 

(N 

in 

m 
m 

■ amplitude of the off-diagonal matrix elements, and a nonzero critical value A c 

o : 

is necessary to magnetize the ground-state. We extend these results by numer- 



ical studies of A for standard tight-binding models which indicate a regime 



■ of intermediate disorder where off-diagonal fluctuations should play an im- 

a: 

portant role for ground-state magnetization. We also emphasize the presence 
of strong correlations between minimal eigenvalues of different magnetization 
5_i ■ blocks which further reduce the probability of a nonzero ground-state spin. 

PACS numbers: 73.23.-b, 7I.I0.-w, 7I.24.+q, 75.10.Lp 
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Itinerant ferromagnetism results from the interplay between the Pauli principle and the 
electronic interactions. The exchange energy can be minimized when the fermionic antisym- 
metry is supported by the spatial wavefunction and when the resulting energy gain exceeds 
the corresponding loss in kinetic energy, this implies spin alignment and a nonzero ground- 
state magnetization. Because of the locality of the Pauli principle, ferromagnetism is often 
studied in the framework of the Hubbard model for which a simple Hartree-Fock picture pre- 
dicts ferromagnetism to occur when the Stoner criterion || is satisfied U c p{ep) = 1, where 
p(€f) is the one-particle density of states at the Fermi energy. Within perturbation theory, 
the effective interaction strength is enhanced by the presence of disorder || and this had led 
to recent predictions of the occurence of a ferromagnetic instability already below the Stoner 
threshold ||. However as we recently pointed out, there is a competing effect of the interac- 
tion manifested through the fluctuations of off-diagonal interaction matrix elements which 
suppresses magnetization [|IJ. Indeed as is known from studies of models in nuclear physics, 
strong off-diagonal fluctuations can determine the variance of the Many-Body Density of 
States (MBDOS) |J. When one includes the spin degree of freedom within such models, 
the Hamiltonian splits into blocks of different total magnetization o~ z and one finds that 
these fluctuations result in a broader MBDOS for the block with smallest magnetization, 
increasing thereby the probability of finding a non magnetized ground-state; therefore the 
average exchange effect must compete not just with the kinetic energy but with off-diagonal 
fluctuations. 

To shortly summarize the results presented in 0, we start from a Hamiltonian for n spin- 
1/2 particles distributed over m/2 one-particle orbitals e a G [— m/2;m/2] (the mean level 
spacing is then fixed as A = 1 with spin degeneracy) 

This Hamiltonian is a generic model of interacting fermions expressed in the basis 
of Slater determinants constructed from the eigenstates ip a of the corresponding free 
fermions model Hq. The interaction matrix elements are given by U2'n — I drdr'U(r — 
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r')^l l (r)ip*p(r')'ip 1 (r)'il)5{r') , where U(f — r') is the interaction potential. Disorder results in 
a random character of the ipa's, leading to fluctuations in U^'t around their average value. 
Only diagonal matrix elements have a nonzero average leading to mean-field charge-charge 
and spin-spin diagonal interactions. We neglect in (1) the mean-field charge-charge contri- 
bution since it has no influence on the magnetization and assume gaussian-distributed U^'p 
: P(U2'b) e~^ Ua <^ ^ 2U . s a = J2st °L s <*s,tc a ,t are spin operators and the parameter A > 



is given by A = (U%p)/ (2RMS [C^JJj, i.e. it is half the ratio of the average exchange 
interaction to the root mean square value of the off-diagonal matrix elements. Finally, the 
interaction commutes with the total magnetization so that the second sum on the right-hand 
side of (1) is restricted to spin indices satisfying s + s' = t + 1' . The Hamiltonian acquires a 
block structure where blocks are labelled by their total magnetization. Each block's size is 
given in term of binomial coefficients as N(a z ) = f^f-o-z) Q+o-J' Most importantly, the 
two-body interaction we consider connects only Slater determinants which differ by at most 
two one-particle indices and each magnetization block's connectivity (the number of nonzero 
matrix elements per row) K(o~ z ) is a monotonously decreasing function of o~ z . 
Similarly as in Q we estimate the variance of each magnetization block's MBDOS as 

tt^ E BlA" - « K (°')& ( 2 ) 
where Hj j = (I\H\J) and \I) is a Slater determinant. Hence each block's bandwidth 
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goes as y K(a z )U. Moreover, earlier work tells us that in such models as (1), the MBDOS is 
approximately gaussian with vanishing corrections in the dilute limit 1 < n C m [|. We can 
therefore conclude that the full MBDOS is a sum of approximately gaussian contributions 
from each spin block with a variance proportional to the corresponding connectivity. The 
latter is a monotonously decreasing function of o~ z , hence the broadest MBDOS corresponds 
to the minimally magnetized block and the ground state will be found in this block with 
increased probability. For approximately gaussian distributions one may assume that the 
tail of each block's MBDOS scales with its variance (up to a numerical factor (3). Neglecting 
contributions from if , the typical spin gap can be estimated (for A = 0) as 



A" w (3U[^K{\a mm \) - ^K(\a min \ + 1)] (3) 

Switching on the mean-field spin-spin interaction Hs induces shifts of the block's MBDOS 
by energies of at most —\Ua 2 z and the spin gap becomes on average 

(A.) = A? - XU (4) 

Results presented in M indicate only a very weak m-dependence of . It has however a 
strong even-odd dependence : A^ pa AU (8U) for even (odd) number n of electrons. This 
effect is nevertheless compensated by a twice larger antiferromagnetic shift between the 
two lowest magnetized blocks for odd n (a z = 1/2 and 3/2) than for even n [o z = and 
1). Hence we introduced A = (3 — (— l) n )A/2 in (§) which results in a parity independent 
behavior of A. Numerically we found a critical value of A c ~ 2 to have a non-negligible 
probability of magnetizing the ground-state. A larger value ~ 4 is then necessary to have 
(A s ) = resulting in a 50% probability of nonvanishing magnetization. 
These points have been studied in details in [|l| . Here we first address the fact that equation 
(|) gives the average value of the spin gap, based on the assumption that the lowest levels 
Eq and Ei in the two blocks with lowest magnetization (they define the spin gap A s = E\ — 
Eq) have uncorrelated distributions. To check this we study the following cross-correlation 
function 

= m- mm -my (5) 

It is easily seen that if E and E\ fluctuate independently around their disorder averages (E ) 
and (-Ei), R{E , E{) — > whereas strong correlations between E — (E ) and E\ — {E\) - i.e. if 
a large (small) value of one of them corresponds to a large (small) value of the other - implies 
R(Eq,Ei) pa 1. Fig.l shows the [/-dependence of R(Eq,Ei) for the Hamiltonian (1) and 
7i=4, 5 and 6 electrons on m = 10 orbitals (the function is independent on A). The numerical 
data show indeed strong correlations between extremal levels of the two lowest magnetized 
blocks. Even though R(E Q , Ex) is somehow reduced by Hjj, correlations still remain quite 
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strong R(E ,Ei) > 0.8 up to very large U = 10. This is understandable if we remember 
that all blocks, having K(cr z )N(a z ) nonzero matrix elements, are constructed out of the 
same set of only (m(m — l)/2) 2 different two-body interaction matrix elements. Extremal 
eigenvalues in the distribution are then due to special realizations of the latter inside the 
blocks. These realizations are presumably not very different in blocks with consecutive 
magnetization which results in strong eigenvalues correlations. The existence of these strong 
correlations is confirmed by the gap distribution P(A^) for Hu plotted on Fig. 2. Indeed 
the probability to find a zero or negative gap is vanishingly small (J^_ 00 P(A^)d(A^) = 
for n = 5 and < 0.002 for n = 4 and 6), whereas even at large average gap, it would not 
vanish in the case of two independent distributions for E and E\. Fig. 2 indicates an 
effective repulsion between extremal eigenvalues of different random matrices, built however 
from the same (small) set of matrix elements. This repulsion results in a further decrease of 
the probability for a magnetized ground-state. Note that the even-odd dependence of the 
average of these distributions is removed due to the horizontal scaling A u s j{A u s ). 
We now turn our attention to the microscopic computation of the magnetization parameter A 
for realistic solid-state models. Specifically, we concentrate on the two-dimensional Anderson 
lattice 

H = VJ2 c{ s c hS + ]T W4, s ^s (6) 
<i;i> » 

where restricts the sum to nearest neighbors, Wj G [— W/2; W/2] and W is the disorder 
strength. Also we consider different interaction potentials of the form 

U(f- ?) = U 5(r- ?) + Ui/ \r - r'| (7) 

Fig. 3 shows the disorder dependence of A, for the Hubbard interaction case Ui — and 
different linear system size L=10, 20, 50 and 80. The data have been obtained from averages 
over 30 wavefunctions in the middle of the band E = and for 10 (L = 80) to 200 (L = 10) 
disorder realizations. Clearly, three regimes are distinguishable. (/) At low disorder, the 
one-electron dynamics undergoes a crossover from ballistic to diffusive regime as the linear 



system size is increased beyond the elastic mean- free path l e ~ 50(V/W) 2 . In the ballistic 
regime l e 3> L, wavefunctions are plane-waves. In this case, a Hubbard interaction gives 



reached, one expects A ~ A/(A/g) ~ g (g is the conductance) |||7|]. In the regime of 
intermediate disorder, both off-diagonal fluctuations and exchange are increased by disorder, 
and apparently they compensate each other, resulting in a L-independent A ~ 4, i.e. close 
to the critical value to have a vanishing average gap Therefore we can expect that 
in this regime, and even for very strong interaction, the probability for finite ground-state 
magnetization does not exceed 50%. Note that for W/V = 5, the localization length is 
slightly below the largest system size considered and one can expect finite-size effects to 
play only a minor role. In the regime of strong disorder, one-particle wavefunctions 

are strongly localized on fewer and fewer sites, so that many of the interaction matrix 
elements are much smaller than U, the off-diagonal fluctuations are sharply reduced and 
again exchange dominates. Note that eventually, the latter disappears also, but at a lower 
rate than the fluctuations. 

We finally evaluate the influence of the interaction range by considering Ui ^ 0. The 
average exchange interaction term is explicitly given by 



(■■•)a :/ a is an average taken over one-particle wavefunctions close to the Fermi level. Due to 
their orthogonality, taking this average over the full set of wavefunctions gives a ^-function so 
that only the on-site term contributes. This last procedure is exact if the one-body dynamics 
is described by Random Matrix Theory (RMT) for which the structure of the eigenstates is 
homogeneous all through the spectrum. In this case - which is relevant for irregular quantum 
dots with chaotic scattering at the boundary || - only the on-site part of the interaction 
contributes to the exchange, whereas it is expectable that increasing the strength U\ of 
the long-range terms leads to stronger off-diagonal fluctuations, and a reduction of A. For 
the Anderson model however, one-particle wave-functions are different from RMT so that 



A ~ L 2 , since the RMS [u^j 



L- 4 and (U%%) ~ L 



2 , whereas once the diffusive regime is 




(8) 
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the average over wavefunctions close to the Fermi level leads only to a more or less sharply 
peaked function of (r — r') : there are also contributions to the exchange from the long-range 
terms, but still we expect that the average damps them with respect to their contribution 
to off-diagonal fluctuations (This damping of course depends on the disorder strength.) so 
that A should decrease with increasing interaction range. The validity of this reasoning is 
illustrated on Fig. 4 where we plot the evolution of A for different disorders as the long-range 
part of the interaction becomes more and more important. Clearly, A is increased by the 
screening of the electron-electron interaction, and therefore the Hubbard results presented 
on Fig. 3 give an upper boundary for A. One thus expects the demagnetizing effect described 
here and in [|]] to be more efficient at low filling when the screening length exceeds the elastic 
mean free path. 
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FIG. 1. Cross-correlation function (5) between the lowest energy levels in the two blocks with 
lowest magnetization for 1000 realisations of the Hamiltonian (1) with n=4 (circles), 5 (squares) 
and 6 (diamonds) particles on m = 10 orbitals. 
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monds) particles on m = 10 orbitals. The distributions have been constructed from 20000 realisa- 
tions of the Hamiltonian 
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FIG. 3. Magnetization parameter A vs. disorder strength W/V, for a Hubbard interaction and 
a two-dimensional 10x10 (full circles), 20x20 (empty squares), 50x50 (full diamonds) and 80x80 
(empty triangles) Anderson lattice. One clearly differentiates three regimes : (I) At small disorder, 
A increases due to a crossover from ballistic to diffusive behaviour (see text). (II) At intermediate 
disorder, exchange and fluctuations compensate each other so that A is size-independent. (Ill) At 
large disorder one-body states are strongly localized over very few sites, which kills the off-diagonal 
fluctuations faster than the exchange and the latter dominates again. 
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FIG. 4. Effect of the interaction range on the magnetization parameter A for a two-dimensional 
15x15 Anderson lattice and ratio Uq/U\ = 1 (full circles), 4 (empty squares), 9 (full diamonds) 
and oo (Hubbard interaction - empty triangles). As expected, an increase in the interaction range 
leads to a stronger increase of the fluctuations than of the exchange, resulting in a lowering of A. 
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